Spin-polarized tunneling as a probe of (Ga,Mn)As electronic properties 
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We present magnetic and tunnel transport properties of (Ga,Mn)As/(In,Ga) As/(Ga,Mn) As struc- 
ture before and after adequate annealing procedure. The conjugate increase of magnetization and 
tunnel magnetoresistance obtained after annealing is shown to be associated to the increase of both 
exchange energy A exc h and hole concentration by reduction of the Mn interstitial atom in the 
top magnetic electrode. Through a 6x6 band k.p model, we established general phase diagrams 
of tunneling magnetoresistance (TMR) and tunneling anisotropic magnetoresistance (TAMR) vs. 
(Ga,Mn)As Fermi energy (Ef) and spin-splitting parameter (Bg). This allows to give a rough 
estimation of the exchange energy A e:cc / l =6BG — 120 meV and hole concentration p~1.10 20 cm -3 of 
(Ga,Mn)As and beyond gives the general trend of TMR and TAMR vs. the selected hole band 
involved in the tunneling transport. 

PACS numbers: 72.25.Dc ; 75.47.-m ; 75.50.Pp 



I. INTRODUCTION 

In the field of spintronics, the p-type ferromagnetic 
semiconductor (Ga,Mn)As offers many advantages to 
study tunnel magnetotransport properties when used 
as an electrode. The complexity of the transport 
mechanisms associated with spin-orbit coupled states 
make this material a powerful means for finding novel 
effects and provides new challenges for theoretical 
understandings. This includes tunnel magnetoresistance 
(TMR) across single and double barriers^ tunnel 
anisotropic magnetoresistance (TAMR), 3,4 Coulomb 
blockade anisotropic magnetoresistance^ and current 
induced magnetization switching. 6,7 However one of the 
main limitation of this p type material for spintronic 
integration is the relatively low Curie temperature. 
Through low temperature annealing trcatcmcnt after 
growth, Curie temperatures of 173 K can be obtained. 8 
Elimination of interstitial manganese atoms which are 
double donors and which couple antiferromagnetically 
with the manganese atom in substitutionnal position, is 
mainly invoked. These atoms diffuse towards the surface 
to form either a MnO-^-or a MnN 11 layer, depending 
on annealing conditions. 

In this paper we describe the effect of anneal- 
ing on the magnetic and electric properties of a 
(Ga,Mn)As/(In,Ga)As/(Ga,Mn)As tunnel junction. We 
have focused our report on this single structure even 
though other junctions with different Mn concentrations 
and ferromagnetic layer thicknesses were studied, lead- 
ing to the same general conclusion^ In the first part we 
detail the effect of annealing on magnetization measure- 
ments and confirm observations made on a (Ga,Mn)As 
trilayer structure with a GaAs barrierji 3 - The second 
part presents the results obtained on junctions fabri- 
cated by optical lithography and describes the behaviour 
of Resistance Area (R.A) product, Tunnel Magnetoresis- 



tance (TMR) and Tunnel Anisotropic Magnetoresistance 
(TAMR) through annealing. In the last part, a general 
interpretation of the data behaviour from both magnetic 
and electric measurements, is given through a 6x6 band 
k.p model of the tunneling transport. Two important 
parameters are identified, the Fermi energy and the spin 
splitting parameter Be introduced in the framework of 
the Zener model in the mean-field approach^ 



II. EXPERIMENTAL RESULTS 

Ga .926Mn .o74As (80nm)/ Ino.25Gao.75As (6nm)/ 
Gao.926Mno.o74As (15nm) structure is grown by molec- 
ular beam epitaxy at 250 °C on a p-doped GaAs buffer 
layer (p=2-10 19 cm~ 3 ). Annealing treatment has been 
realized at 250 °C in a nitrogen atmosphere during 
1 hour. The annealing was performed on a whole 
piece of 5x5 mm 2 for magnetic measurements whereas 
realised on patterned junctions for electrical experiments. 

Figure 1 presents magnetization behaviour before and 
after annealing by SQUID (Superconducting Quantum 
Interference Device) measurements. The two step mag- 
netization reversal along [100] axis at 10 K is due to the 
consecutive reversal of the two magnetic layers [Fig. QJi] . 
As a function of annealing, three important characteris- 
tics, consistent with the reduction of Mn interstitials in 
the top magnetic layer, can be extracted from those mea- 
surements : a large decrease of the coercivity He as well 
as an increase of the magnetic moment M5 and of the 
Curie temperature Tc- Concerning the variation of the 
He, magnetization study pointed out that this decrease 
may be due to an elimination of interstitial manganese 
(double donors) pinning center^ The resulting increase 
of the carrier concentration may also contribute to the 
decrease of the anisotropy field. 14 Considering that only 
the top layer is affected by annealing, through its linear 
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FIG. 1: (Color online) (a) Magnetization measurements vs. 
magnetic field at 10 K before and after annealing along [100] 
direction; (b) Magnetization measurements as a function of 
the temperature before and after annealing in a field of 500 
Oe. 

dependence on the magnetization saturation value, 
the spin splitting parameter increases from 17 meV 
(before annealing) to 24 meV (after annealing). The 
values of the spin splitting were estimated through the 
relationship Bq= a ^^ s derived from the mean field 
theory, where is the Fermi Liquid parameter and f3 
the p-d exchange integral 

The observed Curie temperature are in good agreement 
with those found on thicker magnetic layers confirming 
that layer width larger than 50 nm should still have a 
high concentration of manganese interstitialsi 10 i 16 In 
the present case, the Curie temperature goes from 55 
K to 122 K [Fig. QJb)]. Due to the higher magnetic 
moment of the top layer, the behaviour of the thin 
bottom layer is hidden, supporting the results that only 
the top layer properties change. A further confirmation 
that annealing does not act on the bottom layer comes 
from Auger measurements, not presented here. A strong 
manganese accumulation at the top of the surface is 
measured, whereas no obvious change in the bottom 
layer is observed, already put forward by Chiba et 
alJ£ Capping (Ga,Mn)As layer by a simple GaAs layer 
which width exceeds 5 nm, does not improve the Curie 
temperature of the simple magnetic layer— and does 
though support our results. The formation of a p-n 
junction avoiding the migration of interstitial n-type 
manganese has been suggested. 18 

Magnetic tunnel junctions have been patterned by 
optical lithography (size of the junctions were between 
8 and 128 ^m 2 ). With standard dc technique the 
resistance of the junctions is measured at 3K and at 
low bias (1 mV) in the CPP (Current Perpendicular to 
Plane) regime. Non-linear I(V) curve indicates that a 6 
nm (In,Ga)As layer still acts as a barrierii The reason 
is on the one hand that the Mn acceptor level in GaAs 
leads to a positive band offset in (Ga,Mn)As compared 
to GaAs. On the other hand the well-known As antisites 
incorporated during the low growth process should 
probably govern the pinning of the Fermi level and a 
higher barrier than the simple Mn acceptor state may 
be expected. In figure HJa) we note an increase of TMR 
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FIG. 2: (Color online) (a) Tunnel magnetoresistance mea- 
surements as a function of the magnetic field at 1 mV and 3K 
for a 128 fim 2 junction, (b) Tunnel magnetoresistance mea- 
surements as a function of Resistance. Area product at 3 K for 
4 (un)annealed junctions, (c) Tunnel magnetoresistance at 1 
mV as a function of the temperature before and after anneal- 
ing, (d) Tunnel anisotropic magnetoresistance measurements 
as a function of the magnetic field at 1 mV and 3K. 



from 30% before annealing to 120% after annealing on 
a 128 /jm 2 junction (along [100] direction) while R.A 
product decreases from 0.047 to 0.003 tt. cm 2 . Lowering 
R.A. product should be related to a change in the Fermi 
energy which involve a reduction of the barrier height 
or the barrier width and then must be associated to 
an increase of the hole concentration. In addition, as 
already observed on magnetic properties, the coercive 
field of the top magnetic layer changes after annealing; 
the difference of those values between magnetic and 
transport measurements is related to size effects. The 
same behaviour has been observed on all 4 measured 
junctions [Fig. [2jb)]: Whereas TMR values lay between 
30% and 90% before annealing, an homogenization of 
the values occurs after annealing where TMR ranges 
between 110% and 130%. No assymetry of TMR 
between positive and negative applied bias has been 
measured after annealing. 

However, we must emphasize that magnetic properties 
derive from the whole magnetic layers (volume effect), 
whereas electric properties should mainly depend on the 
interfaces between the tunnel barrier and the electrodes. 
It results that evaluating the change of the electronic 
properties for each magnetic layers from transport 
measurements appears more complex than in the case of 
magnetic experiments. Nevertheless, some conclusions 
can be drawn from TMR measurements vs. temperature, 
taking into account that the elementary process is a 
spin-conservative direct tunneling, i.e. the evolution of 
TMR with temperature is directly linked to the effective 
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carrier spin polarization of the ferromagnetic layer. 19 
We note that the effective temperature at which TMR 
cancels remains unchanged after annealing [Fig[2jc)]. 
This feature comes from the magnetic properties of the 
bottom electrode which are not modified after annealing 
(T c ~ 55 K). The drop of TMR around 15 K before and 
after annealing is related to the quick variation of the 
coercive field of the thin magnetic layer as a function of 
the temperature^ 

On the other hand, how behaves the tunnel anisotropic 
magnetoresistance (TAMR)? TAMR generally traduces a 
variation of resistance vs. the crystalline orientation of 
the electrode magnetization. In this case, this originates 
from the anisotropy of the valence band of (Ga,Mn)As. 
Careful attention was paid on the resistance difference 
when the magnetization is aligned along [100] (in plane 
magnetization) or [001] (out of plane magnetization) 
which leads to maximum TAMR effect in our samples. 
In a saturating field of 6 kOe variations are almost equal 
to 10-15% before and after annealing [Fig. H(d)], in good 
agreement with experiments obtained on a ZnSe barrier. - 
When driving experiments in the plane of the layer resis- 
tance variations as small as 4% were recorded. 

Combination of magnetization and transport measure- 
ments let us therefore presume that the spin splitting 
and the Fermi energy play an important role in tunnel- 
ing transport. The influence of those parameters will be 
discussed now through a 6x6 band k.p modelisation of 
spin-orbit coupled states tunneling transport. 

III. THEORETICAL MODEL 

Our calculations of the transmission coefficient are 
based on the multiband transfer matrix technique devel- 
oped in details by Pethukov et al£^, Brey et al£2- and 
Krstajic et al~ and applied to the hole 6x6 valence 
band k.p Hamiltonian H^. Added to the Kohn-Luttinger 
kinetic Hamiltonian, this includes a p-d exchange term 
introduced by the interaction between the localized Mn 
magnetization and the holes derived in the mean-field 
approximation thus giving: 

H h = -(71 + 4 72 )fc 2 + 672 J2 L ^ k l + 

a 

+673 X! ( L <* L P + L (j L a )k a k fj + \ so ~tl) + 6B G fnS (1) 

equivalent to the one proposed by Dietl et al^ and Abol- 
farth et al^. Here, a — {x,y,z}, L a are I — 1 angular 
momentum operators, is the vectorial spin operator, 
m the unit magnetization vector and 7$ are Luttinger 
parameters of the host semiconductor GaAs. 6Bg rep- 
resents the spin-splitting between the heavy holes at the 
Tg point like originally introduced by Dietl et al^ We 
do not take explicitly into account the stress hamiltonian 
which is shown to give the same qualitive conclusions. 



o — GaAs = 1110.25030.75^ 
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FIG. 3: (Color online) Calculated Resistance. Area (RA) 
product of the trilayer structure as a function of the band 
offset ds between the ferromagnetic semiconductor and a 6 
nm barrier of GaAs or (In,Ga)As. Insets : (Bottom) RA. 
product as a function of the (In,Ga)As barrier width d. (Top) 
Valence band profile of the considered heterostructure. 



To derive the transmission coefficient, the boundary 
conditions to match at each interface are^i: 

i) the continuity of the 6 components of the envelope 
function according to V>rt+Err r «,«Vv= E„' *n,n'V$ 
where the subscript t n n > (r ni n) refer to the respective 
transmission (reflection) amplitude from incident (n) , re- 
flected (n) and transmitted (n ) waves together with 

ii) the continuity of the 6 components of the cur- 
rent wavevector according to J4 , n+J2n Tn ,nJ''^ii = 
Sn' nl Jv> where, in the k.p approach, the current 
operator in the z direction writes J = j- . 

Concerning the heterostructure itself, the valence band 
offset (VBO), d B , between (Ga,Mn)As and (In,Ga)As 
fixes the effective barrier height <p according to &b=- 
"Ep+4> where E^~-0,18 eV is the Fermi level within 
(Ga,Mn)As calculated from the top of the (Ga,Mn)As 
valence band [Inset Fig[3]. On figure [31 we present the 
calculated R.A product vs. the respective valence band 
offset using standard Landauer formula of conductance 
for 6 nm GaAs and Ino.25Gao.75As barriers. Although the 
VBO between Gao.926Mn .o74As and Ino.25Gao.75As is 
still unknown, recent photoemission spectra determined 
the barrier height 4> between (Ga,Mn)As and GaAs to 
450 meVj2£ in agreement with our k.p model considering 
a R.A product approaching ~ 10 -3 il.cm 2 [Fig. [5] and 
like obtained experimentally by Chiba et The rela- 
tive small band offset between valence band of GaAs and 
(In,Ga)As inferior to 50 meV)2& makes then such value 
of </>~ 450 meV a plausible order of magnitude for the 
effective barrier height for Ino.25Gao.75As matching with 
the R.A product after annealing. However, in the present 
case, the real value of 4> may vary depending on: 

i) the nature and density of the dangling bonds at 
the interfaces promoted by the low-temperature growth 
procedure^ 
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ii) the local density and position in energy of ionized 
defects such as As antisites in the barrier which strongly 
influences the valence band bending of the whole het- 
erostructure. 

However, surprisingly, we have noticed that TMR re- 
mains quasi insensitive to the barrier height </> (not 
shown). Let us then focus on the phase diagrams TMR 
(Ef , B G ) and TAMR (E F , B G ) established from the pre- 
ceeding model of k.p tunnel conduction. Figure 4 displays 
both TMR and TAMR vs. the spin splitting parameter 
B G and the Fermi energy Ep of the ferromagnetic semi- 
conductor whose properties are assumed to be identical 
at both sides of the (In,Ga)As barrier. The zero energy 
for Ef corresponds here to the top of the valence band for 
(Ga,Mn)As in its paramagnetic phase (Bg—0). Also are 
plotted on phase diagrams 3 different lines correspond- 
ing to constant carrier concentrations of 1-10 20 cm -3 , 
3.5T0 20 cm~ 3 and 5T0 20 cm -3 , as well as the energy of 
the four first bands at the center of the Brillouin zone. 
Let us first emphasize on the general trends for TMR and 
TAMR from such diagrams. 



IV. DISCUSSION 

High TMR values, up to several hundred percents, 
can be expected either for spin splitting values larger 
than several tens of meV or for low carrier concentra- 
tion, that is when only the first subband is involved in 
the tunnelling transport. This corresponds to a quasi half 
metallic character for (Ga,Mn)As. Starting from the first 
subband and increasing the carrier concentration to fill 
the consecutive lower subbands (n=2,3,4), up and down 
spin populations start to mix up, leading to a decrease of 
TMR. For high carrier concentration (n=4), small TMR 
is expected which may anticipate difficulties to concili- 
ate high Curie temperature and large TMR effects. We 
specify that for low values of spin splitting and Fermi 
energy, ferromagnetic phase induced by carrier dereal- 
ization may not exist (top right corner of the diagram) 
which is not taken into account in our k.p modelisation 
(propagative envelope wave function). In the same man- 
ner, we cannot reproduce metal- insulator transition in 
the tunneling transport, responsible for the large TAMR 
obtained in in-plane geometry^^ 

What about TAMR signal ? We can firstly note a 
possible change of sign for TAMR on crossing the third 
subband. The first subband clearly gives a negative 
contribution to TAMR. This originates from the pre- 
dominant heavy hole character of such band, an in-plane 
magnetization allowing, through off diagonal compo- 
nents, a possible heavy to light hole conversion, and 
then a larger transmission through the barrier^ This 
argument is reversed for the second and third subbands 
with the results that TAMR becomes positive when 
n=2 and n=3 subbands are dominant in the tunneling 
transport. We can point out that a change of TAMR 
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FIG. 4: (Color online) Tunnel magnetoresistance values (a) 
and tunnel anisotropic magnetoresistance values (b) repre- 
sented as a function of the Fermi and spin splitting energy 
for a 6 nm (In,Ga)As barrier with a band offset of 450 meV. 
White lines represent the 4 bands at the center of the Bril- 
louin zone. Gray lines indicate the Fermi energy for different 
hole concentrations. 



sign was already observed on a Zener-Esaki diode 3 ^ as 
well as theoretically established through tight-binding 
treatment^ Reducing the hole concentration through 
hydrogenation technique should give the possibility to 
probe this possible crossover from positive to negative 
TAMR.^ 2 - 

Concerning our own experiments, taking into account 
conjugate TMR and TAMR values obtained before and 
after annealing, one can roughly evaluate the projection 
of the corresponding signals trajectories in the \Ep, B G ] 
plane followed during annealing [FigHJ- A good qualita- 
tive agreement can be found even though symmetrical 
junctions were simulated in order to restrict the number 
of parameters. 

Evaluating directly the interfacial spin splitting 
from the mean field theory appears difficult since the 
interfacial magnetic properties are hardly accessible. 
However, when using the estimated B G value of the 
top magnetic electrode (before and after annealing) 
a good qualitative agreement can be found for TMR 
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and TAMR, as illustrated by the trajectory in figure [4] 
between point 1 (before annealing) and point 2 (after 
annealing). A more refine calculation including two 
different Bq after annealing should be required to draw 
definite quantitative conclusion. 

We are now going to discuss the hole concentration 
derived from these diagrams. TMR and TAMR values 
obtained before annealing are well reproduced for a hole 
concentration approaching 10 20 cm~ 3 , in good agreement 
with the one measured for single (Ga,Mn)As layer and 
already reported^ The annealing procedure has for 
effect to i) remove Mn interstitial atoms, ii) increasing 
carrier concentration and Hi) reduce the effective barrier 
height even if the valence band position is expected 
to rise due to an increase of the average exchange 
energy (Be). The large reduction of the R.A product 
together with the increase of TMR are consistent with 
such assumption. Nevertheless, the hole concentration 
extracted after annealing from the phase diagram 
~l,7.10 20 cm -3 appear to be weak compared to the one 
reported in the literature and derived from Hall effect 
measurements. The existence of a possible concentration 
gradient can be at the origin of such discrepency. Also 
can be invoked, a reduction of the hole concentration 
at the interfaces with the barrier due to a significant 
charge transfer between p-type (Ga,Mn)As and n-type 
(In,Ga)As (excess of As antisites)i 28 ' 34 



properties of the top magnetic layer, ensuring an in- 
crease of the effective magnetization and a significant 
enhancement of the tunnel magnetoresistance. The 
confrontation between experiments and modelisation 
performed within a 6x6 band k.p treatment vs. intrinsic 
(Ga,Mn)As parameters (hole filling, exchange energy) 
allowed a rough estimation of the average exchange 
interactions and carrier concentration in (Ga,Mn)As at 
the interface with the barrier. We point out that while 
the magnitude of TMR appears very sensitive to both 
parameters (Bq and E^), the TAMR variation is limited 
to several tens of percent but may change sign crossing 
from upper to lower (Ga,Mn)As subbands. As a final 
conclusion, we think that this reduced parameter model 
gives a good qualitative agreement of the tunneling 
transport and enables to extract the fundamentals of 
TMR and TAMR processes involving tunnel transport of 
spin-orbit couple state. In order to go further and draw 
more quantitative information, a perfect control and 
knowledge of the carrier density seems to be necessary. 
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V. CONCLUSION 

In summary we have shown that annealing a 
(Ga,Mn)As-based tunnel junction mainly affects the 
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